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Abstract 

The physics of a row of topphng dominoes is discussed. In particular the 
forces between the fahing dominoes are analyzed and with this knowledge, the 
effect of friction has been incorporated. A set of limiting situations is discussed 
in detail, such as the limit of thin dominoes, which allows a full and explicit 
analytical solution. The propagation speed of the domino effect is calculated 
for various spatial separations. Also a formula is given, which gives explicitly 
the main dependence of the speed as function of the domino width, height and 
interspacing. 

1 Introduction 

Patterns formed by toppling dominoes are not only a spectacular view, but their dy- 
namics is also a nice illustration of the mechanics of solid bodies. One can study the 
problem on different levels. Walker P] gives a qualitative discussion. Banks j2] con- 
siders the row of toppling dominoes as a sequence of independent events: one domino 
undergoes a free fall, till it hits the next one, which then falls independently of the 
others, and so on. He assumes that in the collision the linear momentum along the sup- 
porting table is transmitted. This is a naive viewpoint, but it has the advantage that 
the calculation can be carried out analytically. A much more thorough treatment has 
been given by D. E. Shaw [Sj. His aim is to show that the problem is a nice illustration 
of computer aided instruction in mechanics. He introduces the basic feature that the 
domino, after having struck the next one, keeps pushing on it. So the collision is com- 
pletely inelastic. In this way a train develops of dominoes leaning on each other and 
pushing the head of the train. One may see this as an elementary demonstration of a 
propagating soliton, separating the fallen dominoes from the still upright ones. Indeed 
Shaw's treatment is a simple example how to handle holonomous constraints in a com- 
puter program describing the soliton. As collision law he takes conservation of angular 
momentum. We will demonstrate, by analyzing the forces between the dominoes, that 
this is not accurate. The correction has a substantial influence on the solition speed, 
even more important than the inclusion of friction, which becomes possible when the 
forces between the dominoes are known. 

The setting is a long row of identical and perfect dominoes of height h, thickness 
d and interspacing s. In order to make the problem tractable we assume that the 
dominoes only rotate (and e.g. do not slip on the supporting table). Their fall is due to 



the gravitational force, with acceleration g. The combination provides a velocity 
scale and it comes as a multiplicative factor in the soliton speed. Typical parameters of 
the problem are the aspect ratio d/h, which is determined by the type of dominoes used, 
and the ratio s/h, which can be easily varied in an experiment. Another characteristic 
of the dominoes is their mutual friction coefficient which is a small number (~ 0.2). 
The first domino gets a gentle push, such that it topples and makes a "free rotation" 
till it strikes the second. After the collision the two fall together till they struck the 
third and so forth. So we get a succession of rotations and collisions, the two processes 
being governed by different dynamical laws. Without friction the rotation conserves 
energy, while the constraints exclude the energy to be conserved in the collision. In 
fact this is the main dissipative element, more than the inclusion of friction. 

The goal is to find the dependence of the soliton speed on the interdistance s/h. 
In the beginning this speed depends on the initial push, but after a while a stationary 
pattern develops: a propagating soliton with upright dominoes in front and toppled 
dominoes behind. The determination of the forces between the dominoes requires that 
we first briefiy outline the analysis of Shaw. Then we analyze the forces between the 
dominoes. Knowing these we make the collision law more precise. With the proper 
rotation and collision laws we give the equations for the fully developed solitons. The 
next point is the introduction of friction and the calculation of its effect on the soliton 
speed. As illustration we discuss the limit of thin dominoes (i — > 0, with permits for 
small interseparations a complete analytical solution. Finally we present our results 
for the asymptotic soliton speed for various values of the friction and compare them 
with some experiments. We also give an explicit formula, which displays the main 
dependence of the soliton speed on the parameters of the problem. The paper closes 
with a discussion of the results and the assumptions that we have made. 



2 Constraints on the Motion 

The basic observation is that domino i pushes over domino z + 1 and remains in contact 
afterwards. So after the contact of i with i + 1 the motion of i is constrained by the 
motion of i + 1. Therefore we can take the tilt angle 9n of the foremost falling domino, 
as the only independent mechanical variable (see Fig. . Simple goniometry tells that 

hsm{6i — Oi^i) = {s + d) cos6'j+i — d. (1) 

To see this relation it helps to displace domino i + 1 parallel to itself, till its bottom line 
points at the rotation axis of domino i (see Fig. Q). By this relation one can express 
the tilt angle 9i in terms of the next 6'j+i and so on, such that all preceding tilt angles 
are expressed in terms of 9^. The recursion defines 9i as a function of 6'„ of the form 

9i=pn-^{9n): (2) 

i.e. the functional dependence on the angle of the head of the train depends only on 
the distance n — i. The functions Pj{9) satisfy 

= + arcsin J(£±i?)£25pW^) , (3) 

with the starting function po{9) = 9. They are defined on the interval < 9 < 9c, 
where 9c is the angle of rotation at which the head of the train hits the next domino 

9c = axcsin{s/h). (4) 
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We will call 6c the angular distance. From the picture it is clear that the functions 
are bounded by the value 6^0, which is the angle for which the right hand side of 
vanishes ^ 

cos 6100 = — —:■ (5) 
s + a 

600 is the angle at which the dominoes are stacked against each other at the end of the 
train. We call 6*00 the stacking angle. 




Figure 1: Successive dominoes. The tilt angle 9i is taken with respect to the vertical. 
In the rectangular traingle ABC the top angle is a = 6'j — 6*^+1, the hypotenuse has 
the length h and the base BC the length (s + d) cos 6i^i — d. Expressing this base in 
the hypotenuse and the top angle yields relation (^. In the picture the tilt angle of 
the head of chain On has reached its final value 9c = aj:csm{s/h). The first domino has 
almost reached the stacking angle 6^0 ■ The normal force fi and the friction force /i/j 
that domino i exerts on i + 1 are also indicated. 



The picture shows that the functions pj {6) are monotonically increasing functions. 
They become flatter and flatter with the index j and converge to the value 600 (at 
least not too close to the maximum separation s = h, see Section [TUI). The functions 
are strongly interrelated, not only by the defining recursion The angle 6i can be 
calculated from the head of the train 6n by Pn-i but also from an arbitrary intermediate 
6k by Pk^i- This implies 

Pn-i{0) =Pk-t{Pn-k{0)), e.g. Pj{0) =Pj-i{pi{0)). (6) 

One easily sees that pi{0) = 9c. Therefore one has 

pM=Pj-iipM)=Pj-M), (7) 

a property that will be used later on several times. 

An inmediate consequence of (0) is the expression for the angular velocities Ui = 
d9/dt in terms of a;„. From the chain rule of differentiation we find 

dOi dOji . , 
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with 

..(.)^^. (9) 

Computationally it is easier to calculate the Wj recursively. Differentiation of © with 
respect to 9^ yields 

^jW^u,,-.w(i- ^ '^+''"""^"'. 1. (10) 

\ h COS [pj{6) - pj_i{^^^ ' 



Another useful relation follows from differentiation of the second relation ^ 

wj{9) = w,^Me))wiie) w,io) = w,^iiec), (ii) 

since pi{0) = 9c and wi(0) = 1. 

3 Rotation Equations 

Without friction, the motion between two collisions is governed by conservation of 
energy, which consists out of a potential and a kinetic part. The potential part derives 
from the combined height of the center of mass of the falling dominoes, for which we 
take the dimensionless quantity 

n 

HniOn) = EK^^ + (d/h) sin ft,]. (12) 

i 

The kinetic part is given by the rotational energy, for which holds 

n 

K„ien,cOn) = {I/2)Y^ul I=il/3)mih' + £), (13) 

i 

where / is the angular moment of inertia with respect to the rotation axis and m is 
the mass of the dominoes. We write the total energy as 



En = ^mghcn = ^mgh (^Hn{On) + ~Jl^^(^^) ' (1^) 



where the dimensionless effective moment of inertia I{6n) is defined as 

n 

In{en) = E ^l(^n)- (15) 

j 

We have factored out mgh/2 in (fT^ as it is an irrelevant energy scale. This has the 
advantage that the expression between brackets is dimensionless. The factor I/mgh 

I _h{i + dyh')_ ^^g^ 



mgh 3g 

provides a time scale that can be incorporated in u;„. From now on we put this factor 
equal to unity in the formulae and remember its value when we convert dimensionless 
velocities to real velocities. 
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We see (II 4p as the defining expression for ujn as function of 9^ 



As mentioned e„ is a constant during interval n. So we can solve the temporal behavior 
of 9n from the equation 

= UJn[9n). (18) 

The initial value for 9n is and the final value equals the rotational distance 9c- The 
duration of the time interval where n is the head of the chain, follows by integration 

Jo UJn[9n) 

In this time interval the soliton has advanced a distance s + d. The ratio (s + d)/tn 
gives the soliton speed, when the head of the train is at n. In order to integrate 
the equations of motion ()18|) we must have a value for e„ which basically amounts to 
finding an initial value co'n(O) as one sees from ()14|) . In the next section we outline how 
to calculate successively the oJniS^). 

Putting all ingredients together we obtain the asymptotic sohton speed fas as 

^ ( %^ lim 1. (20) 

In this formula the time t„ is computed from the dimensionless equations (setting 
I/mgh equal to 1). 



4 The Collision Law, first version 

We now investigate what happens when domino n hits n + 1. In a very short time 
domino n + 1 accumulates an angular velocity u;„+i(0). The change in cOn+i takes place 
while the tilt angles of the falling dominoes hardly change. Shaw [Sj postulates that 
the total angular momentum of the system is unchanged during the collision. This is 
not self-evident and we comment on it in Section IHl Before the collision we have the 
angular momentum 

n 

L^ = Y.w,{9c)un{9c). (21) 

i 

After the collision we have 

n+l 

Ln+l = Y.^M^n+M- (22) 

i 

Equating these two expressions yields the relation 

n n+l 

o;„+i(0) = ujn{9c) Y.^j{9c)/ ^M- (23) 

j j 

With the aid of this value we compute the total energy e„+i and the next integration 
can be started. For the first time interval holds cq = 1 + t^o(O) since only the zeroth 
domino is involved and it starts in upright position with angular velocity uJo{0). The 
value of ujq{0) has no influence on the asymptotic behavior. After a sufficient number 
of time intervals, a stationary soliton develops. 
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5 Forces between the Dominoes 



Conservation of energy requires the dominoes to slide frictionless over each other. 
Before we can introduce friction we have to take a closer look at the forces between the 
falling dominoes. Without friction the force which i exerts on i + 1 is perpendicular to 
the surface of i + 1 with a magnitude fi (see Fig. H]). Consider to begin with the head 
of the train n. Domino n feels the gravitational pull with a torque T„ 

= (sin^„-(rf//i)cos0„)/2, (24) 

and a torque from domino n — 1 equal to the force times the moment arm with 
respect to the rotation point of n. The equation of motion for n becomes 

^ = T„ + fn-i h[cos{9n-i - On) - {s + d) slu 6*^+1 ]. (25) 

Domino n — 1 feels, beside the gravitational pull T„_i, a torque from n which slows it 
down and a torque from n — 2 which speeds it up. Generally the equation for domino 
i has the form ^ 

-J^ = Ti + fi-itti-i - fibi. (26) 

The coefficients of the torques follow from the geometry shown in Fig. ^ 

= hcos{9i — Oi^i) — {s + d) sin 6*4+1, = hcos{9i — 9i+i). (27) 

Note that the first equation is just a special case with /„ = 0. Another interesting 
features is that < 6j. So i gains less from i — 1 than i — 1 looses to i. Therefore 
dominoes, falling concertedly, gain less angular momentum than if they would fall 
independently. This will have a consequence on the application of conservation of 
angular momentum in the collision process. We come back on this issue in the next 
section. 

We can eliminate the forces from the equation by multiplying with vq = 1 and 
the general equation with r„_j and chosing the values of rj such that 



(ro = 1), or r„_, = \[^. (28) 



"n-j j=i '-'j 

Then adding all the equations gives 



duji 
It 



Y^[fi^irn-iai_i- fiVn-i-iai] =0. (29) 



Now observe that the recursion for the Tj is identical to that of the Wj as given in ()10j). 
With ro = 1 we may identify Vj = Wj. It means that if we multiply ()29|) with Un and 
replace Tn-iUJn by Ui, we recover the conservation of energy in the form 

|iE-?-E-.r.. (30) 

It is not difficult to write the sum of the torques as the derivative with respect to time 
of the potential energy, thereby casting the conservation of energy in the standard 
form. So if conservation of energy holds, the elimination of the forces is superfluous. 
However, equation (j29p is more general and we use it in the treatment of friction. 
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6 The Collision, second version 



We have assumed that in the coUision of the head of chain n with the next domino n + 1 
conserves angular momentum. Having a more detailed picture of forces between the 
sliding dominoes we reconsider this assumption. In this section without friction and 
in Section |H1 with friction. The idea is that in the collision domino n, exerts a impulse 
on n + 1 and vice versa with opposite sign. In other words: one has to integrate the 
equations of motion of the previous section over such a short time that the positions 
do not change, but that the velocities accumulate a finite difference. However, not 
only the jump in velocity propagates downwards, also the impulses have to propagate 
downwards in order to realize these jumps. Denoting the impulses by capital F's, 
domino i receives Fi from i + 1 and from i — 1. So we get for the jumps in the 
rotational velocity 



'^n.+l(0) — Fnttn, 
Wi{0) a;„+l(0) - Woi9c) i^n{Oc) = K-l ^n-l - K bn, 

. w„,+i_i(0) t^„+i(0) - w„_i(6'c) uJniOc) = Fi_i ai_i - Fibi. 



(31) 



The functions and bi are the same as those defined in (j27|) . If we would have = bi 
we could add all equations and indeed find that the angular total angular momentum 
is conserved in the collision. But only a„ = bn since 6n+i = 0. The impulse Fi can be 
eliminated in the same way as before by multiplying the ith equation with r„_|_i_j and 
adding them up. For the coefficient of ujn+i{0) we get 

n+1 n+1 

r„+i_iw;„+i_i(0) = J2 '^i^i(O) = -Jn+i, (32) 

i j=0 

and for the coefficient of ujn{0c) one finds with (fTUj) 

n n n 

Yrn+i-iWn^i{d) = ^r„+i_i w„,+i_i(0) = ^^rjiUjiO) = J„+i - 1- (33) 

i i j=l 

As general relation we get 

Jn+l C^n+l(0) = ( Jn+l - 1) t^n(6'c)- (34) 

In our frictionless case rj = wj and therefore Jn+i = /„+i(0). So the desired relation 
reads 

/„+l(0) CUn+m = (/n+l(0) - 1) UJnie,) = ^{9,) UJniO,). (35) 

We have added the last equality since it smells as a conservation of angular momentum 
using the effective angular moment of inertia I{0). This inertia moment is however 
linked to the energy and not to the angular momentum. The true angular momentum 
conservation is given in Section EJ It is also not conservation of kinetic energy. Then 
the squares of the angular velocities would have to enter. The difference with the 
earher relation is that the sum involves the squares of the w's. This has a notable 
influence on the asymptotic velocity. 
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7 Fully Developed Solitons 



After a sufficient number of rotations and collisions a stationary state sets in. Then 
we may identify in the collision law the entry uJn+i{0) with ^^.^(O). This allows to solve 
for the stationary C(j„(0). We use (jllll to relate the effective moments of inertia 

n n—1 

ini9c) = E^K^c) = T.^'+M + w'M) = W)-wlio) + wl{e,). (36) 

j=0 j=Q 

For large n the last term vanishes and we may drop the n dependence in /„. So 

1(6,) = /(O) - 1. (37) 

The collision laws thus may be asymptotically written as, 

/(0)cc;„(0) = [/(0)-l]^„(ej. (38) 

The rotation is governed by the conservation of energy, which we write as 

m ulie) + H„,{e) = /(O) u^O) + H^{0). (39) 

We can use Q to relate the height function Hn{9c) to its value at 6^ = 0. 

" d d 

Hn{Oc) = Y.[^ospj{e,) + - sinpjie,)] = Hn{0) - 1 + cospnidc) + ^ sinp„(^e)- (40) 

The limiting value of p„ is the stacking angle 6oo Therefore the difference between the 
initial and the final potential energy reads 

H{0)-H{e,) = 1 - cose^ - ^sin^oo = P{h,d,s). (41) 

We have introduced the function P as the loss in potential energy in the soliton motion. 
It is the difference between an upright domino and a stacked domino at angle 6*00- The 
functional form reads explicitly 

It is clear that the domino effect does not exist if P is negative, because a domino 
tilted at the stacking angle has a higher potential energy than an upright domino. 
We use (1411) in the conservation law for the energy, taken at 6 = 6c 

1(6^) iuUec) - /(O) iuUO) = P{K d, s). (43) 
Solving a;„(0) and a;„(6'c) from (jHHj) and (jlSj) yields 

u;l{0) = P{h, d, s) ^M^, ujHOc) = P{h, d, s) j^^. (44) 



By and large yP sets the scale for the rotation velocity. The dependence on J(0) is 
rather weak. For large /(O) it drops out. The minimum value of /(O) is 2 which is 
reached for large separations. 
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8 Friction 

After all this groundwork it is relatively simple to introduce friction. Let us start with 
the equation of motion (j26|) . Friction adds a force parallel to the surface of z + 1. For 
the strength of the friction force we assume the law of Amonton-Coulomb jl] 

/friction = fJ'f, (45) 

where / is the corresponding perpendicular force. Inclusion of friction means that the 
coefficients and bi pick up a frictional component. The associated torques follow 
from the geometry of Fig. ^ So the values of the Oj and bi change to 

tti = h cos{9i — 9i+i) — {s + d) sin 6*4+1 — f^d, 

(46) 

bi = hcos{9i - 9i+i) + fihsm{9i - 9i+i). 

Then we may eliminate the forces as before, which again leads to (|29|). But we cannot 
identify any longer rj with Wi. In order to use (j29|) we must express the accelerations 
dui/dt in the head of chain dun/dt. This follows from differentiating (jH)) 

Wn-.{9n)^ + Vr..0) Ul (47) 



y^^e) = (4J 



dt -""-'^""^ dt 

with Vi given by 

dwj{9) 

lt9n 

The Vj can be calculated from the recursion relation, that follows from differentiating 
(fTn|) . Clearly the recursion starts with f o = (see (jEj)). 
Next we insert ()47|) into ()26|) and obtain 




The equation can be transformed into a differential equation for duJn/d9n by dividing 
dnOl) by u;„ = d9n/dt 

We use this equation to find as function of 6'„ and then (fTH|) again to calculate the 
duration of the time between two collisions. 

The inclusion of friction in the collision law is even simpler, since relation 
remains valid, but now with the definitions pUj) for and bi. 



9 Thin Dominoes 

Sometimes limits help to understand the general behaviour. One of the parameters, 
which has played sofar a modest role, is the aspect ratio d/h. In our formulae it is 
perfectly possible to take this ratio 0. In practice infinitely thin dominoes are a bit 
weird, because with paperthin dominoes one has e.g. to worry about friction with the 
air. In this limit we can vary s/h over the full range from to 1. In Fig. |21we have 
plotted the asymptotic velocity as function of the separation s/h. The curve is rather 
fiat with a gradual drop-off towards the large separtions. We discuss here the two 
limits where the separation goes to and where it approaches its maximum s = h. 
Both offer some insight in the overall behavior. 
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Figure 2: The asymptotic soliton velocity as function of the separation s/h in the thin 
domino hmit. Also is plotted the result of Banks in the same limit. 



9.1 Infinitesimal Separation 

If the dominoes are narrowly separated, the head of chain rotates only over a small 
angle 9c = arcsin(s//i) c:^ s/h and the collisions will rapidly succeed each other. The 
number of dominoes with a tilt angle 6i between and 7i/2 becomes very large and 
slowly varying with the index i. So a continuum description is appropriate. We first 
focus on the dependence of 6i{6n) on the index i and later comment on the dependence 
on the weak variation with 6n (which is confined to the small interval < 6n < Oc)- 
We take as coordinate x the distance of domino i. 

X = i s/h (51) 

and use v = ns/h for the position of the head of the train. Then 

e^ = e{x), e^+i = 6{x + dx), (52) 

with dx = s/h. So for d = Q and s/h ^ the constraint (P) becomes 

sin[ 6{x) — 9{x + dx) ] = dx cos 9{x + rfx), (53) 
leading to the differential equation 

de{x) 



dx 

which has the solution 



-cos6'(x), (54) 



sin6'(x) = tanh(z/ — x) or 9{x) = arcsin(tanh(z/ — x)). (55) 
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Here we have used the boundary condition that 0(z/) = 0. Not surprisingly we find 
that the shape of the tih angles is a function of the difference with respect to the head 
of the train. The above expression gives the shape of the soliton. 

Next we comment on the dependence of this profile on the angle On- As mentioned 
it can be only weak as the interval for On is narrow. Thus it suffices to know a few 
derivatives and for that, the interpretation (jlU|) is useful. The behavior of Wj in the 
continuum limit, follows from the differential form of the recursion relation 



dw{x) 
dx 



sin 0{x) w{x) 



with the solution 



1 



wix] 



cosO{x). 



(56) 



(57) 



cosh(i/ — x) 

Note that, not unexpectedly, the form of w{x) follows also from that of 0{x) by differ- 
entiation with respect to u. Similarly the expression for Vj, as given by (jlHj) . can be 
obtained from differentiation of (jFfj) with respect to u 

tanh(z/ — x) dw{x) 



vix] 



cosh(z/ 



x] 



dx 



(5J 



What still is needed is the propagation velocity of the soliton, or in the present 
language: how fast n or u moves with time. As the foremost domino rotates over 
a small angle Oc — s/h, the head of train covers the distance s/h with the rotation 
velocity Un- So the propagation speed equals Un- As before, Un has to be distilled 
from the laws of rotation and collision. Since this section is mainly for illustration, we 
restrict ourselves to the frictionless case. 

In the collision law (|35p we encounter uJn{Oc) and u;„_|_i(0). Both are linked to 
a;„(0) = uj{u) by 



(jJn{Oc) = + ^^s/h, 

OUn 



a;„+i(0) = uj{iy) + ^^/h. 



(59) 



For the derivative with respect to On-, we can take advantage of the form ()50p which 
directly gives this derivative. We use that rj = Wi in the frictionless case. The sums 
can be performed explicitly in the continuum limit using (|57p and (|58|) 
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S " 



1 ^ 



dx 



X] 



dx 



cosh 

tanh(z/ — x) 
2 cosh(z/ — x) 



dxw(x) — — 
dx 



Therefore the equation for dun/dOn becomes 



tanh V, 



cosh V — 1 
2 cosh V 



tanh V 



duJn cosh u — 1 1 1 2 

tanh u — — = — h - tanh u Un- 

oOn 2 cosh U Ur, 2 



With (|59p the collision equation has the form 



dcu 



du 



(60) 



(61) 



(62) 
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Using (jUUI) we get, to first order in s/h, 



tann I/— — = ti;(z/) + tanh z/— . (63) 

OO^ Oh' 

Next we substitute (jHTjl and we obtain the following differential equation for ^^(z/) 

dujiv) , . A ,9 X cosh I/ — 1 1 , 

tanhz/— ^ = cj(z/)(-tanh2i/ - 1) + — . (64) 

dv ^ '^2 ' 2 cosh z/ uj{v) ^ ' 

This awful looking differential equation has a simple solution 

cj^(z/) = "^^^ „ — (cosh V — 1 — log cosh v) . (65) 
sinh V 

We have chosen the integration constant such that oj{v) vanishes for v = 0. It starts 
as 

tu(z/) ~ Z//2V2 + ---, (66) 



and it saturates exponentially fast to the value c<j(oo) = 1 (leading to Vas = y/3gh). 
Thus we have obtained in the continuum limit a full and explicit solution. It may serve 
as an illustration for the general discrete case. 



9.2 Maximal Separation 

On the other side, near maximal separation s ^ h, also a simplification occurs. Here 
the number of dominoes involved in the train is restricted to a few. The head of the 
train rotates over almost 7r/2 before it strikes the next domino. So one comes close 
to the picture of Banks [2j in which the toppling of the dominoes is a succession of 
independent events. There is however a difference resulting from the constraint (0). 
Inmediate after the collision, the dominoes n and n + 1 rotate with equal velocity 
<-^n{Oc) = i^n+i(0). This is a consequence of the fact that after the collision, one still 
has 6n+i = 0. Thus we find for the energy after the collision 

e„+i = l + 2^2^i(0). (67) 

All other dominoes have fallen down and domino n + 1 is still upright (the 1 in (jUTjl ). 
Once n+1 starts rotating, the value of Un rapidly drops down to 0. Inspecting recursion 
firU]) . with po{dn+i) = ^n+i and Pi{On+i) = On, one sees that the factor cos(^^„ — 9n+i) — 
cos(^c) is very close to 0. The ratio approaches 

s sin(6'„+i) _^ s sin(6'n+i) _ s_ 
h cos{9n — On+i) h cos(7r/2 — 9n+l) h 

So Wi — s> and indeed domino n comes to a halt; is has to, since it has reached the fioor. 
This has an effect on the moment of inertia I{6n+i) defined in (fT3j) . Inmediately after 
the collision the value of I{9n+i) equals 2, being the sum of w^+i = 1 and 10^ = 1. A 
small angle further it has dropped to 1, since Wn drops to 0. As the energy is conserved 
the kinetic energy of domino n is transferred to domino n + 1. So Un+i rises by a factor 
\/2 in a short interval. Therefore we start the integration of the time after this sudden 
increase, using the conservation law for the energy 

ul^^iOn+i) + cos^„+i = 1 + 2a;2^i(0). (69) 
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In particular we have the relation for 9c — ti/2 

^^+i(^c) = l + 2a;^+i(0). (70) 

The collision law for this degenerate case becomes 

o;„+i(0) = a;„(ej/2. (71) 

Note that since the Wj are either 1 or 0, there is no difference between the proposal by 
Shaw (see (1221)) ^'^d ours (see (jHSl))- 

The stationary state is obtained by the identification ujn{9c) = u!n+i{dc). Combining 
flTOj) and (f7T|) then yields ciJ„+i(0) = l/\/2. Thus the time integral for the interval 
becomes in the stationary state 

t = / d9 , = 1.37 (72) 

Jo v2 — cos^ 

The reciprocal yields the asymptotic soliton speed v^s = 0.73 * -\/3 = 1.26. 

The story of thin dominoes gives a warning on the numerical integration scheme. 
For small separations we need many intervals before the asymptotic behaviour has set 
in. On the other hand we do not need many points in the integration for the time of 
a rotation. For wide separations it is the opposite: only a few intervals are needed for 
the asymptotic behavior, but we have to perform the time integration with care. The 
factor I{6) in the energy law is rapidly varying for small 6. So we need many points 
for small 6 to be accurate. 

In Fig. 121 we have also plotted the curve due to Banks [2j in the limit of thin 
dominoes. The difference is due to the collision law, for which Banks takes 

UJn+l{0) = COsdcUJn{0c). (73) 

The factor cos 6c accounts for the horizontal component of the linear momentum. For 
large separations this gives quite a different value, since the transmission of linear 
momentum becomes inefficient. For small separations the transmission is nearly perfect 
(as becomes our collision law). The conservation of energy of a single rotating domino 
reads 

uj^{9) + cose = LU^ {0) + l. (74) 
For the stationary state we insert (f7H|) into (f7i|) and get 

9//, \ 1 ~ cos6'c 
^ (^-) = 1 • (75) 

1 — COS^ dc 

For small 9c, the value to {9c) approaches 1/v^, which is again substantially smaller 
than our limiting value 1. The reason is that the dominoes, which keep leaning onto 
each other and onto the head of the train, speed up the soliton. 



10 Calculations and Limitations 

For the frictionless case we can use the formulae of Section [71 i.e. first calculate the 
asymptotic value of uj{0) and then integrate the rotation equation to find the time 
between two collisions and thus the asymptotic soliton speed. With friction we must 
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s/h 


Shaw 


frictionless 




= 0.1 


/X = 0.2 




= 0.3 


0.1 


3 


64568 


2 


23469 


1 


82095 


1.51695 


1 


28221 


0.2 


3 


06742 


1 


95534 


1 


66019 


1.43423 


1 


25452 


0.3 


2 


74686 


1 


82515 


1 


56987 


1.37279 


1 


21498 


0.4 


2 


50849 


1 


74231 


1 


50459 


1.32193 


1 


17605 


0.5 


2 


30183 


1 


67865 


1 


44771 


1.27272 


1 


13420 


0.6 


2 


10337 


1 


62447 


1 


39204 


1.22009 


1 


08609 


0.7 


1 


89899 


1 


57824 


1 


33347 


1.15958 


1 


02745 


0.8 


1 


68680 


1 


53779 


1 


26196 


1.07970 





94711 


0.9 


1 


47984 


1 


47984 


1 


15267 


0.95822 





82681 



Table 1: The asymptotic sohton speed {d/h = 0.179) for the colhsion law of Shaw and 
for various degrees of friction with the collision law ()34|) . 



iteratively find a;(0), by trying a value of a;(0), then solve equation (jSUj) for u!{9c) and 
finally apply the collision law in order to see whether we come back to our trial a;(0). 
A form of iteration is to start the train with one domino and an arbitrary initial a;o(0) 
and let the train grow longer such that an asymptotic pattern develops. 

In Table ^ we have summarized the results. The thickness to height ratio is set at 
d/h = 0.179 since this is the only value on which experiments [H] are reported. The first 
column gives the separation s/h, the second the soliton speed using Shaw's collision law 
and the third gives the results for ours (j35|) . In the subsequent columns the influence 
of the friction is indicated. Note that the reduction of the speed due to the change of 
the collision law is larger than that of modest friction. The curves corresponding to 
these values are shown in Fig. El which also contains the experiments of Maclachlan 
et al. They suggest that the soliton speed diverges for short distances, while we find a 
maximum. Their values seem to correspond best with the friction coefficient /i = 0.3. 
We found empirically the value /i = 0.2, by estimating the angle of the supporting 
table at which dominoes start to slide over each other. 

In order to make the behavior of fully developed solitons more transparant, we may 
introduce, for frictionless dominoes, the average 

1 1 fOc de 

(76) 



with uo{6) the solution of ()39|). This average is a number close to 1/ \/P (with P defined 
in ()42|) ). since the integrand varies from a value slightly larger than l/\/P to a value 
slightly less than 1 / \/P. Then we get for the asymptotic soliton speed the formula 

g(/z,rf,s)^J^^, (77) 



^ \/Pih,d,s) 
where the factor Q is given by 



Qih, d, s) = ( ^ ] (^+^) ^nKd,s) ^ ^^^^ 

' ' yl + d^/h^j ha,icsm{s/h) 
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Figure 3: The influence of friction on tfie asymptotic soliton speed for the aspect ratio 
d/h = 0.179. The dots are the experimental values of Maclachlan et al. jH]. 



Here we have reinstalled the factor I /mgh in order to include in this formula, all 
the factors that contribute to the velocity. The factor Q is shown as function of s/h 
for various d/h in Fig. |3] One may consider Q as the main factor determining the 
dependence of the soliton speed on the parameters of the problem. The fraction in 
(f77j) is a refinement which requires a detailed calculation. We found that this fraction 
is virtually independent of the aspect ratio d/h. It stays close to 1 for the the major 
part of the range of practical separations. Only around the already "unworkable" 
separation s/h = 0.9 the value has increased some 10%. A good indicator for the 
behavior is the curve for the frictionless thin dominoes which is the product of the 
fraction and Q = \/3 s / {h aTcsm{s / h) . 

We mentioned in Section [3 that the function P as given by (j42|l has to be positive 
for the existence of the domino effect. This gives a bound on the minimal distance s/h, 
which can be cast in the form 

£ > ^Wi>f 

h i-{d/hf ^ ' 

Separations smaller than the value of (f?^ do not show the domino effect and slightly 
above that limit the train has difficulty to develop. The reason is that after a while, 
too many dominoes of the train get tilt angles, which have a higher potential energy 
than an upright domino. Ultimately the fraction of these dominoes in the train looses 
out against the dominoes at the end of the train, which are tilted at the stacking angle 
(with a potential energy lower than an upright domino). One can overcome this barrier 
by starting with an unreasonble high initial uJo{0). So ff^i is the true theoretical limit, 
but in practice the domino effect will not start for slightly larger values of s/h. 

Another limitation of the theory is at the other side. The dominoes at the end of 
the train are tilted at the stacking angle 6*00 provided the height h is sufficiently large. 
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Figure 4: The function Q{s, d, h) as defined in (f7S|) for various d/h. 



The condition is 

/i' > (s + df - d\ (80) 
For smaller h the dominoes fall flat on the supporting table. (IHU]) is satisfied for 



s/h < ^1 + {d/hy - d/h. 



Beyond this value the train is actually shorter than blind application of the formulae 
would suggest. It is not so interesting to sort out what precisely happens if (jHT|) is 
violated, since then the no-slip condition for the dominoes is highly questionable. For 
such wide separations the force on the struck domino has hardly a torque to rotate 
it. It rather induces the rotation axis to slide along the table. In fact, as a practical 
limitation, we look to the height of impact. If it is above the center of mass of the 
struck domino, it will start to rotate and below that value, it may slip if the friction 
with the supporting table is not large enough. This criterion yields the limit to the 
distance 

s/h <Vs/2 = 0.87, (82) 

which is already a large separation, not far from the limit set by ()81|) for d/h = 0.179. 

For very thin dominoes is hardly a limitation. However, fl79|l and ()82|1 form 
a window of separations for the existence of the domino effect, which depends on the 
thickness d/h. This window narrows down to zero and the domino effect disappears 
for 

h^ < hd^ + Ad^/V3 , or d/h < 0.3787 (83) 

This estimate comes close to the one given by Freericks j5j. Friction also makes the 
excluded interval larger. For fi = 0.2 we have not found a domino effect for s/h < 0.07, 
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which is, for d/h = 0.179, about 7 times the theoretical hmit. So our estimate (j83|) for 
the upper thickness is still too optimistic. 

11 Discussion 

We have studied the toppling of a row of equally spaced dominoes under the assump- 
tions that the dominoes only rotate and that they keep leaning onto each other after a 
collision with the next one. By and large we follow the treatment of Shaw , who intro- 
duced the constraint (^, which synchronizes the motion of train of toppling dominoes. 
By analyzing the mutual forces between the dominoes, we have corrected his collision 
law and we could also account for the effect of friction between the dominoes. The 
correction of the collision law is more important than the influence of friction, given the 
small friction coefficient between dominoes. The limit of thin dominoes d/h leads 
to a completely tractable model. For large separations we encounter a situation which 
resembles the viewpoint of Banks P] , seeing the toppling as a succession of independent 
events. However his collision law differs substantially from ours and cannot be recon- 
ciled with the force picture that we develop. We give a formula (f77|). which displays 
explicitly the main dependence of the soliton speed on the parameters of the problem. 
The maximum speed which can be reached, appears close to the closest separation for 
which the domino effect exists. 

The assumptions, on which our calculation are based, are the no-slip condition and 
the constraint (^. One can help the no-slip condition by increasing the friction with 
the supporting table (putting them on sandpaper as Walker P does). If the no-slip 
condition is violated, it is the end of the domino effect as the dominoes are kicked over 
with the wrong rotation. We argue that this will happen when the falling domino hits 
the next one below its center of mass. 

The constraint (QJ) is implied by the assumption that the collision is fully inelastic. 
This assumption is supported by slow motion pictures of the effect, which show that 
the dominoes indeed lean onto each other while falling. It is an interesting question 
what happens, if the collision would be less inelastic. The extreme opposite, fully 
elastic collisions, yields an ever increasing soliton speed. A falling domino increases its 
rotation velocity as soon as its center of mass goes down. If this is fully transmitted 
to the next domino, the rotation velocity keeps increasing. In that case friction can 
not play a role since the dominoes do not touch each other. In the less extreme case of 
partially inelastic collisions, the dominoes also rotate without contact, but friction can 
play a role during the collision. As Fig. [T] indicates, friction always rotates the mutual 
impulse such, that the torque on the next one decreases and that the reaction torque 
increases. This will slow down the train and a stationary state can develop. Therefore 
it would be interesting to experiment with dominoes of different making (e.g. steel) to 
see the increase in the soliton speed. 

Acknowledgement. This study was motivated by a question from the Dutch 
National Science Quiz of 2003 (www.nwo.nl/quiz). The author is indebted to Carlo 
Beenakker for drawing his attention to the problem, for supplying the relevant literature 
and for stimulating discussions. 

References 

[1] Jearl Walker, Scientific American, August 1984. 



17 



[2] Robert B. Banks, Towing Icebergs, Falling Dominoes and other adventures in 
Applied Mechanics, Princeton University Press, 1998. 

[3] D. E. Shaw, Am. J. Phys. 46 (1978) 640. 

[4] See e.g. D. Tabor, ASME Journal of Lubrication Technology 103 (1981) 169. 



[5] J. K. Freericks, http: / /www.physics.georgetown.edu~jkf/ class_mech/demol.ps 



[6] B. G. MacLachlan, G. Beaupre, A. B. Cox and L. Gore, Falling Dominoes, SIAM 
Review 25 (1983) 403. 



18 



